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Abstract 

The singlet coupling to the topological charge density in the instanton vacuum, causes 
the instantons and antiinstantons to be screened over distances of the order of 1/2 fm. 
Dilute instanton systems behave as a free gas, while dense instanton systems behave as 
a plasma. The free gas behaviour is favored by a density of 1 fm~ 4 . Owing to the Higgs 
mechanism, the ff mass is heavy (1100 MeV). The vacuum topological susceptibility is 
small (0.07 fm~ 4 ) and consistent with the QCD Ward identity. In the chiral limit, the 
singlet screening vanishes, leading to a four dimensional plasma state with a temperature 
given by the nonet decay constant / (90 MeV). The phase is Debye screened. In the 
presence of matter, the screening is quenched, and a Kosterlitz-Thouless transition may 
occur signaling the restoration of the Z7a(1) symmetry. We suggest to use the spatial 
asymptotics of the static topological charge correlator together with the topological sus- 
ceptibility to probe the interplay between the Ua(1) restoration and the chiral restoration 
in finite temperature QCD. The relevance of these results to the bulk thermodynamics in 
the instanton vacuum is discussed. 
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1. Introduction 

The instanton approach to the QCD ground state has received continuous attention 
in the past few years |2|, || |J. Overall, it provides a fair description of the bulk 
vacuum properties and hadronic correlations. Recent (quenched) lattice simulations || 
have suggested that the bulk structure and hadronic correlations are mildly affected by 
cooling, an indication that semiclassical physics may be operative in the long wavelength 
limit. 

Extensive calculations using instantons have been carried out in the quenched approx- 
imation (equivalently large N c limit), both analytically |], and numerically ||. While 
the numerical calculations were found to be consistent with the cooled lattice simula- 
tions, these results suggest that the instantons in the vacuum state are in a dilute gas 
approximation. To the annoying exception of the ff mass, which comes out to be too 
heavy in the numerical instanton calculations, the overall part of the hadronic spectrum 
is reproduced by fitting the pion decay constant and the fermion condensate. The model 
does not confine, yet it breaks spontaneously chiral symmetry with the generation of a 
momentum dependent constituent mass. 

To understand some of the features of the model, it suffices to recall few thinghs. 
The pseudoscalar nonet follows from symmetry provided that the pion decay constant 
and the fermion condensate are reasonable, given suitable current quark masses. In the 
dilute gas approximation, the vector and isobar correlators are not affected by instantons 
at small distances. The formers because of the self-duality character of the instanton 
configuration, and the latters because of the flavor character of the induced instanton 
interaction. Interesting effects in the nucleon channel have been reported suggesting 
diquark correlations at intermediate distances ||. Since the model lacks confinement, the 
large distance behaviour of these correlators is dominated by constituent quark thresholds. 

A major unknown in instanton calculations resides in the character of the instanton- 
antiinstanton interaction. Earlier calculations based on variational estimates with a sum 
ansatz in the quenched approximation, have led to a reasonable description of the bulk 
structure ||. However, these estimates were found to be quantitativaly sensitive to the 
short distance character of the interaction. The latter was ansatz dependent. More elab- 
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orated numerical analyses, based on stream-line configurations as the 'best' instanton- 
antiinstanton prototypes, have shed question marks on the use of a semiclassical ap- 
proximation altogether ||. However, these configurations hold only in the quenched 
approximation, where the instanton-antiinstanton interaction is dipole-like at large dis- 
tances. They may not be important in the unquenched approximation, where screening 
by charges in the fundamental representation is in effect. 

Semiclassical physics in QCD has been challenged by Greensite |T(J using a com- 
bination of arguments based on large N c , factorization and the master field equation. 
Greensite's objections to instantons (and generically to semiclassical models based on 
statistical ensembles) may be evaded by noting that the large N c limit does not commute 
with the thermodynamical limit. 

In this paper we would like to come back to the annoying problem in the ff channel. 
Since instantons have been professed to cure the C/a(1) problem pTT|] , this discrepancy 



cannot be ignored. Also, in the unquenched approximation, the instanton-antinstanton 
interaction is not dipole-like at large distances, but screened over distance scales of the 
order of 1/2 fm. This screening is caused by the mixing between the singlet and the 
topological density in the vacuum, and is at the origin of the large rf mass, thus resolving 
the Ua(1) problem. In fact this construction was briefly presented in [0 and more recently 



discussed in the context of the strong CP problem Ul2| . We will also show that the 
resulting vacuum topological susceptibility is consistent with the QCD Ward identity. 
In the unquenched approximation, dilute systems of instantons behave as a free gas. 
Throughout, the short distance part of the instanton-antiinstanton interaction will be 
assumed. In the presence of matter the screening decreases. Our approach suggests a 
dual description for the bulk pressure : either as a free meson gas with a heavy rj\ or 
as a free meson gas with a light singlet plus a screened Coulomb gas of instantons and 
antiinstantons. At temperatures for which the screening length becomes comparable to 
the interparticle distance, a Kosterlitz-Thouless transition [[HJ may occur. This transition 
may be related to the Ua(X) restoration in QCD. We suggest to use the large spatial 
asymptotics of the topological charge correlator to account for the restoration of ^a(I) 
in the QCD ground state. We explicitly show that the bulk thermodynamics at low 
temperature is meson dominated, and a rapid crossover in the topological susceptibity 



is expected with the vanishing of the fermion condensate. Throughout, the octet-singlet 
mixing will be ignored and the vacuum angle will be set to zero. 

2. Model 

Few years ago, Nowak, Verbaarschot and I suggested that in the long wavelength 
approximation the salient features of the instanton model to the QCD ground state fol- 
lows from symmetries and anomalies. If we were to denote by n the local densities of 
instantons and antiinstantons (thought about here as quasiparticles), then the partition 
function in the presence of a pseudoscalar source term S a , reads 

Z[S] = fdn* J dK e -w\K\-w l{ K,n]+J^{SK) (1) 
where W[-ftT] is the pseudoscalar effective action 

W[K] = J dS {-{d^f - \^J^ Tr(min) (2) 
and Wi[K, n] is the singlet coupling to the local topological charge {n + — n~) ~ E a ■ B a , 



Wj[K, n] = -i I d A x^-^(n + - n~) K°(x) (3) 
-' J 

We have defined the nonet by K = J2l=oK a X a , with the normalization Tr(A a A & ) = 5 ab . 
The mass matrix m = diag (m, m, m s ) ~ diag (5, 5, 120) MeV. The imaginary part follows 
from the fact that (|3D is T-odd in Minkowski space. The singlet coupling to the topological 
charge is y/2Np/ f, where / ~ 90 MeV is the bare nonet decay constant, and Np the 
number of flavours. It vanishes in the large N c limit as it should. The fermion condenstate 
< ijjijj >~ —(255 MeV) 3 , for an instanton density < N > /V 4 ~ 1 fm -4 . We arrived at ([!]) 
by integrating over a random system of instantons and antiinstantons and then bosonizing 
the multi-t' Hooft vertices. We note that eventhough Wi is subleading in 1/N C , since the 
instanton-antiinstanton density < N > /V4 is of order N c , a statistical averaging of this 
part leads effects of order N® (see below), thus comparable to the free meson action W. 
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If we were to drop this term, the present analysis is totally in agreement with the detailed 
analysis of Diakonov and Petrov [Q, in the two flavour case. This term reflects on the 
fermionic induced interactions and will be important for the discussion below. We have 
explicitly disregarded the scalar fields for convenience. The latters are relevant for the 
issue of the scale anomaly and the bulk compressibility. We refer to for a discussion of 
these issues. We stress that the fields K are auxiliary (integrated over), while the physical 
meson sources are represented by S. The auxiliary fields arise from the bosonization of 
the multiflavour quark interactions triggered by a random instanton gas. The a posteriori 
legitimacy of the random gas approximation will be discussed below. 

To leading order in N c and in the long wavelength limit, the fermionic correlation 
functions follow from (^). Specifically 



(^ 75 A a V(x)^75A>(0)) ~ -N c Tr(\ a S+(x)\ a S F (-x 



I < W > 2 ( SHnZjS] \ 

P \5S«(x)5S"(0)) s=0 1 ) 

where Sp(x) is the constituent quark propagator, and the trace is over the spin variables. 
It suffices to know, that asymptotically Sp(x) ~ M 2 e _M ' x '/|x|, where M ~ m + M(0) ~ 
m + 345 MeV 0]. Similar expressions can be derived in other channels. The presence of 
the two-constituent quark cut in (Q) reflects on the lack of confinement in the model. It is 
interesting to compare expression (ffl) with a similar expression derived in two dimensional 



QCD for large N c [14]. In the latter, the two-fermion cut is infrared sensitive and cancels 
precisely against the infrared sensitive one-gluon exchange graph. The result is a sum of 
mesonic poles only. This cancellation makes explicit use of Ward identities in Feynman 
graphs. It is quantum, and thus absent in the semiclassical approach followed here. 

For finite N c the fermions feed back on the instantons and vice-versa. A straightfor- 
ward way to incorporate this, is to use ([[]) for finite iV c . This is an approximation, as 
subleading parts in the multi-'t Hooft interactions were disregarded in (|I]). We expect 
these parts to renormalize the fermion condensate and the nonet decay constant. With 
this in mind and in the long wavelength limit, we can approximate the instantons and 
anti-instantons by point like structures (the effects of their sizes will be discussed below) 
and define the topological charge to be (N = N + + iV_) 



(n + -n )(x) =^2QiSi(x - Xi) (5) 

i=l 

with Q± = ±1. Thus, the interaction term between instantons in the long wavelength 
limit is mostly triggered by the singlet K°, through their topological charge 



W J = ^Z^ y £Q i K (x t ) (6) 

We note that the combination (iy/fK°) plays the role of a four dimensional Coulomb 
potential (the i for Euclidean), with Coulomb charges = Qi^2N F / f. This point will 
be further discussed below. If we denote by z± the fugacity of the quasiparticles, then 
substituting (@) in © gives (N = N + + AL) 

N + N- N / 

Z[ S ] = Y: II / I d \K] e~ W[K]+ ^^ Q * *°<*0+/*»i>(MD (?) 

jv± N +- N -- i=x J J 

The fugacities follows from < N± >= dlnz± Z[Q]/dz±. For a noninteracting system 
lnZ[0] ~ V A (z + + zJ). Thus z± ~ (< N > /2V 4 ). 

In the vacuum (S = 0) we can rewrite (|7|) in two equivalent ways. By summing over 
the instanton-antiinstanton degrees of freedom in the long-wavelength limit, leaving out 
the singlet field. Thus 



Z[0] = J d[K] e - w W + f d4 * 2z -(Vro //) = J d [jq e -s[K] (g) 

which shows that Z[0] is a sum of zero point contributions from a free but light octet of 
pseudoscalars, and a self- interacting but heavy pseudoscalar singlet. This version of Z[0] 
will be relevant for the discussion of matter effects below. Equivalently, we can choose 
to integrate over the singlet K° in (ffi) leaving out the instanton-antiinstanton degrees of 
freedom . Ignoring singlet-octet mixing, the result is 

Z[0] = ( / d'[K]e- w W) £ TFJW1 II / dx i e ~ Nf ~^^ Q,Qi \^ m oi^-^i ) ( 9 ) 
\J J N± N + \NJ. i=l J 
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where m,Q = 2 / '3m 2 K + 1 / 3m 2 is the singlet K° mass, and Ki a Bessel function. The singlet 
K° field is excluded from the mesonic measure, thus the prime. (§) is the product of the 
vacuum partition function of free massive octets and a four dimensional Coulomb gas with 



Yukawa interactions. A similar result was also obtained by Kikuchi and Wudka |L2| . The 
mixing (|7]) causes the instanton-antinstantons in the vacuum to be screened over distances 
of the order of l/m ~ 1/2 fm. Note that the partition function for the Coulomb gas is 
ill-defined Throughout, we will assume that the instantons and antiinstantons 

have a core (smeared). The effects of the smearing will cause the fugacities to renormalize 
by self-interaction terms. In the long wavelength limit, the physics is insensitive to the 
detailed choice of the core. 

The mixing (||) causes the rf excitation to be heavier than the rest of the octet. This 
is just the Higgs mechanism. Indeed, the pseudoscalar correlator correlator follows from 

Zfo/j _ J d[K}e~ W[K]+ 1 22 cos (V2N^K°/f)+ J d^xr,' K°/V2N^ ( 1Q ) 

using @. In the Gaussian approximation (order N® in the action and N c in the correlator) 
the result is 

<5***<«)***<o)> ~ -s.Tr + (0 (^fj) (ID 

where the square of the rj' mass is given by 

m\ =m 2 + 2z -jf (12) 

which is the result quoted in 0. For iVp = 3, / ~ 90 MeV and z ~ 1/2 fm~ 4 we have 
m v t ~ 1100 MeV. The present analysis is consistent with the analysis suggested by Witten 
rS] and Veneziano [fT7| , except for one thingh : z ~ A^ c , so that the contribution of the 



Coulomb gas is z/f 2 ~ N®. The r\-r\ splitting does not vanish in the large iV c limit in the 
instanton vacuum. Should we be concerned ? Yes, if we were to believe that the transition 
from N c = 3 to N c = oo is smooth in QCD (no phase change). Phase changes may occur 
in instanton systems with large values of N c |§. For a dilute system N/V4 ~ 2z ~ 1 fm -4 , 



the interparticle distance is about 1 fm, which is twice the screening length 1/mo ~ 1/2 
fm. In this regime, the instantons and antiinstantons behave as a free noninteracting gas. 
This result shows that the approximations used in || |6], |7| in which the instanton and 
antiinstanton distributions were taken to be totally random is justified in the unquenched 
approximation (finite N c ). Finally, the induced term (|7|) causes the r/' to self-interact. 
The ff quartic interaction term has a strength of order 2z(y/2Np/ f) 4 ~ 10 2 , which is 
strong. We note that this interaction vanishes in the large N c limit. Finally, we note that 
the rf pole in ( |TTD lies above the two-constituent quark cut, which is about 700 MeV, 
leaving the 7/ with a broad width. This problem is generic to models that do not confine. 

3. T = Susceptibilities 

The topological susceptibility of the vacuum can be calculated using 

x (x - y) = QiS 4 (x - Xi) J2 Qjfa(v - Xj)) (13) 

i 3 

since E ■ B ~ (n + — n~). The expectation value (|T3|) follows from the following generating 
functional 

Z [0] = —, r ft / d ^ eQl9(Xl) e " f ^ L - QiQj V^ ) (14) 

by taking d 2 lnZ j d6{x)dd(y) and setting 9 to zero. A straighforward calculation gives 



xO - y) 



P 



(m v , - m ) 5 4 (x - y) - (m, n , - m ) 



2 



'Ki(m^\x - y\Y 



m v > \x 



(15) 



4 F \ 2tt 2 

At large separations, the susceptibility falls off exponentially at a rate given by the phys- 
ical ff mass. This fall off agrees with numerical simulations We note that (plj) is 
analogous to the result in the two dimensional Schwinger model, where the susceptibility 



is given by the correlation function of the electric field [IS 



Xs{x-y) 



Att 2 



< E{x)E{y) >-- 



47T 2 



- y) 



^K (m s \x-y\) 



(16) 



with ms = e/^/n the photon screening mass and Ko a Bessel function. The topological 
susceptibility obeys a zero momentum sum rule. Indeed, Fourier transforming (|15|) gives 
at zero momentum 



2m + m a <ifnf)> ( 2m + m a <if)ifi>\ . . 

Xiq = 0) = = 2ivH 1+ 3 pm% 1 (17) 



in agreement with the QCD Ward identity ]K 



< ijjijj > i nr 



rf 4 x (G(z)G(O)) = ~m^— + I d X <^(^75^(0)> (18) 



where we have used G = (g 2 /16n 2 )E a ■ B a . For the parameters quoted above, x ~ -07 
fm~ 4 , which is small. Note that in the chiral limit, the topological susceptibility (|17D 
vanishes like in the Schwinger model where Xsiq = 0) = 0. 

By analogy with the Schwinger model, the following expectation value 

( e _< Jd 4 <**«G(*)^ =< TT e -«J< J D4 dix *4(«c-aii) > ^ 
i 

can be estimated. Here D4 is a four dimensional volume with D3 = dD^ as a border. 
This expectation value reflects on the amount of topological screening in the vacuum. 
Smilga ]20| has suggested that (fLlp obeys a four- volume law in Yang-Mills theories, and a 



three-surface law in QCD. In the Schwinger model the screening is total, and the analogue 



of (|19|) is the Wilson loop. The latter obeys a perimeter law (Xs(q = 0) = 0). In our case, 
a semiclassical estimate of ( |T9] ) gives 



N r d 

< f[ e~ iQ% d x 64x - Xi) >~ e - 5[K c il (20) 



where is the classical solution to 



-□ + ml)Kl{x) + 2z^^sin {^^K^ = ^={~^ + O d 4 y h(x - y)(21) 
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for a strongly coupled source. While we do not know of a general solution to (21) in 



four dimensions, in the linear approximation (|21|) can be readily solved. In this case, 
the solution is equivalent to taking the second cumulant in (|19D. The latter is just the 
topological susceptibility, 



iV 



< TT e7 iQl ^4 d4x 6 ^ x ~ x ^ >^ e -^x(<?=o) (22) 

i 

Since ( |TTD does not vanish, the second cumulant indicates that ( |H| ) obeys a four-volume 
law. The fall off rate is x/2 ~ 10~ 2 fm -4 . This result is not exact and maybe affected 
by higher cumulants (non-Gaussian fluctuations). We note that the cumulant expansion 
is consistent with the dilute gas approximation. In the chiral limit, the topological sus- 
ceptibility vanishes, and the screening is total. Indeed, for m = m s = fl22|) falls off with 
the border D 3 = <9_D 4 



N 



i 

where m 2 = 4Npz/f 2 with gE the three-normal to the four- volume D^. The fall-off rate 
is given by 

2^ / ( 2 ^ 2 ^ 2 \ ^ 



, )X K rn z v ,--m 2 K --m 2 \ (24) 

which is about 1 fm~ 3 . (|23"D follows from the second cumulant and large three-surfaces 
-D3 m * >> 1. This result is analogous to the result in the Schwinger model. In the chiral 
limit, the instantons and antinstantons are in a strongly coupled plasma-like phase driven 
by a Debye-Huckel equation, with an effective temperature T* ~ / (see below). 

The fluctuations in the number density of instantons and antiinstantons can be ob- 
tained using similar arguments, if we were to note that 

N 

n + + rT = 5 4 (x — %i) (25) 

i 

Since we have ignored the scalars we will only provide a simple estimate in this case. A 
straighforward calculation gives 
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N N 

a(x - y) ~ - a*) 53 &t(y - ^)>c ~ ^z8(x - y) (26) 

i=l i=j 

where only the connected part has been retained. This result shows that the variance in 
the number of particles is Poissonian, and that the compressibility is of order 1 fm -2 . The 
inclusion of the scalars smears the correlations ( ^6|) over the scalar Compton wavelength, 
causing the compressibility to increase. These effects along with the QCD scale Ward 



identities EE will be discussed elsewhere [22 



The present interplay between the instantons and the if with the subsequent screening, 
does not affect the pseudoscalar octet in the long wavelength approximation. Indeed, from 
(H) it follows that the auxiliary meson fields decouple from the topological charge to order 
N®, to the exception of the singlet K° (ignoring singlet-octet mixing). A self-consistent 
calculation would require that instead of using a random gas approximation, we should 



use a screened Coulomb gas approximation, e.g. (14). However, the diluteness of the 
system implies that the effects are small, since the interparticle distance is about twice 
the screening length. Thus we would predict that the fermion condensate and the nonet 
decay constant will not be affected considerably by the screening mechanism. In a way, 
this is good because it implies that most of the calculations performed in the random gas 
approximation reflect correctly on unquenched QCD ||, ^|, |7|]. 



4. T ^ Susceptibilities 

How does the temperature affects the present arguments ? First, consider the correla- 
tion in the topological charge ([13]). At low temperature, the system is meson dominated. 
The temperature enters both through the meson parameters [] as well as the meson dis- 
tributions. Thus (cj n = 2imT) 



+C^(rnl(T) - rn%(T))5 4 (x - y) 



2N f 



(27) 



Tf 2 (T){m 2 v ,{T)-m 2 {T)y_ g ^ n { x °- y °)-J^^)\2-$\ 



4N P 



\x-y\ 



lr The fermion constituent mass becomes temperature dependent |23L \M 
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The static part x( u — 0, x; T) of (27) is dominated by the rf excitation 
X (cu = 0,x;T) = Tj o dx° X (x,0;T) 

_ K(T)-mg(T)) 

2\x\ J 
(28) 

In this form, the QCD Ward identity is fulfilled even at finite temperature, 
f« i n -^x 2m + m,<#> (T) 

J dx X (u J = 0,x,T) = ^— - (29) 

and vanishes at the chiral transition point. This, however, does not necessarily mean 
that the 11^(1) symmetry is restored. The restoration of the latter can be asserted only 
through the large distance behaviour of fl28|). In the Uj\,(l) broken phase, the correlator 
falls off exponentially, while in the symmetric phase it vanishes identically. This point is 
worth probing both on the lattice and in the numerical instanton calculations ||. For 
completeness, we note that the temporal asymptotics of (E8p is free field dominated for 
temperatures T ~ m n > /2tt ~ 150 MeV (0 < a; < 1/T). The compressibility u{T) ~ 
2z(T), and is expected to drop by 1/2 at the chiral transition point, following the drop in 
the gluon condensate [p5H - Debye screening at high temperature causes the electric 
condensate to vanish. In instanton models that would also mean the vanishing of the 
magnetic condensate because of self duality. 

In the Schwinger model, the temperature effects on the topological correlator can be 
calculated exactly. Indeed, the correlation in the topological charge at finite temperature 
is given by 

Xs(x,T) = Lx) - ^1 £ _^_ e ^VV^!l^ (30) 

which is the analogue of (27). The analogue of the static part of the susceptibility (ES|) is 
given by 

Xs (co = 0, x 1 ; T) = jf * X s(x, T) = X ^(x 1 ) - ^e^l^l) (31) 
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whatever T. We observe that eventough integrates to zero, the correlator still falls 
off exponentially. In the Schwinger model mj follows directly from the Ua(1) anomaly 
(bubble graph) and thus is T-independent for all temperatures. This is not the case for 
rrirf as we have discussed above. For completeness, we also note that along the temporal 
direction fl3~0|) is dominated by the free field behaviour for temperatures of the order of 
T ~ ms/2ir (0 < x° < 1/T), 



Xs(x°,T) ~ JL ^ Vl 6( X °) - ^ + g ln( 27 rT|x°|)) (32) 
where V\ is the space length. 



5. Pressure 

At low temperature, the pressure is just given by the one loop effects following from 
(0) and the rest of the mesonic action. Thus 



P + B =P 7r + P^ + P J7 + P 7? / + ... (33) 

where B is the vacuum pressure, and P a the respective mesonic pressures. The dots stand 
for the higher hadron resonances, in agreement with the Gibbs average. The mesons in 
(|33|) carry temperature dependent masses, since the fermion condensate is temperature 



dependent. At low temperature, the effects are small, and the pressures P a in fl33|) are the 
usual black-body contributions. Below the deconfining transition, the vacuum pressure is 
T-independent in a confining theory p7|, to leading order in N c . This is not the case here 
and B receives an unwanted contribution from the free constituent quark loop of order 
N c . 

In writing ( |3"3"D we have ignored the 77' self-interactions, i.e. we have kept the order 
N® terms in the free energy. For finite N c these effects, may not be small. Indeed, strong 
self-interactions may give rise to classical solutions to (0 = iy/fK°) 



—U + m ) <f) = —z \ —f- J ( e / ~ e / I ( 34 ) 
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which shows that <fi (the precursor of the 77') plays the role of a coarse grained Coulomb 
potential in four dimensions. In the chiral limit (mo = 0) (|34]) is the Debye-Huckel 



equation for a four dimensional Coulomb plasma with charges q = yj2Np/f ~ l/v%, 
charge density n ~ zq/Air ~ y/^Nc, and an effective temperature T* ~ / ~ \/W c 0. For large 



N c , the Coulomb plasma is classical, and (|33"D follows. The 77 — rj' mass difference is just 
the Debye screening length. At high temperature, the effects of the dropping masses may 
be important. The derivation, however, becomes less reliable as other effects (quantum 
fluctuations, subleading 1/N C terms, ...) may be important. Note that at high temperature 
the four dimensional Debye-Huckel equation (taken litterally) dimensionally reduces 
to an equation in three dimensions. The effective temperature for the three dimensional 
Coulomb gas is / ~ 0, that is low. The system has a tendency to cluster. This tendency 
is perhaps what has been observed by Ilgenfritz and Shuryak in interacting instanton 
systems with fermion determinants. In a way the high temperature Coulomb gas problem 
discussed here behaves as a low temperature Kosterlitz-Thouless system |13|]. Coulomb 
systems in two dimensions display a dipole phase at low temperature, and a plasma phase 
at high temperature. If such a transition were to occur, it is plausible that the resulting 
phase is Ua(1) symmetric. Since x(0) ~ m < ipip >, it is possible that this phase change 
may coincide with the chiral transition. This point can be elucidated by lattice simulations 
that will calculate both — 0, q = 0,T) and x( u — 0, \x\,T), as suggested above. 

We note that the effects of the instantons in (|33D is implicit, and manifest through the 
larger 77' mass even at finite temperature (ignoring 77' interactions). An equivalent way of 
describing the same physics, is to rewrite ( p3|) as follows 



P + B = P 7r + P^ + P 7? +(P#o+P / ) + ... (35) 

where P^-o is the contribution to the pressure coming from a thermalized system of sin- 
glet mesons with mass m , and Pi the instanton contribution following from the finite 
temperature partition function 

2 A four dimensional plasma, follows from a five dimensional field theory with <\> ~ M 3 / 2 and q ~ M~ x / 2 
where M carries canonical mass dimension. 
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Z/[T ] = ^^ JJ f d * Xie NF 7^)^ Q ^J(^l^m (36) 



f£N + \NJ.£Jn?xi/T 

where q n = (2irnT,q). P r] i = Pro + Pj. The Coulomb gas description referred to above 
is now manifest. This decomposition, allows us to see how matter affects instantons. It 



also shows the artificial character of the decomposition in (35) . At low temperature, the 
condensate (thus mo) and the meson decay constant / do not change appreciably. The 
instanton screening length remains about the same. At high temperature tuq is expected to 
drop, weakening the screening. Our calculations break down when the screening becomes 
comparable to the interparticle distance, i.e. m\ ~< iV > T/V3. Subleading terms in 1/N C 
and the omitted scalars as well as higher lying hadrons are important at high temperature. 



6. Conclusions 

We have given some qualitative arguments indicating the interplay between the 77' 
fluctuations and the instantons in the vacuum state and in matter. At zero temperature, 
the instantons are screened over distances of the order of 1/2 fm, and the 77' acquires a 
mass of the order of 1100 MeV, leading to a natural resolution of the problem, as 



originally suggested by 't Hooft [11]. These results emphasize our earlier work p] and are 



in general agreement with the arguments given by Kikuchi and Wudka fI2 |. The screening 
among the instantons and antiinstantons in the vacuum yields a topological susceptibility 
that is consistent with the QCD Ward identity. For dilute systems of instantons, the 
screening justifies the use of the random gas approximation for finite N c (unquenched 
calculations) [|], ||, [/J. In the instanton model of the vacuum the 77 — 77' mass splitting 
does not vanish in the large N c limit. 

The screening is unaffected by low temperatures, i.e. temperatures for which the pion 
and kaon mass do not change substantially. As expected, the bulk pressure is meson 
dominated. The 77' contribution to the pressure can be viewed as originating either from 
a strongly correlated and massive 77' at temperature T, or from a free light singlet at 
temperature T mixed to a four dimensional Coulomb gas with an effective temperature 
f(T). The latter can be described by a Debye-Huckel equation in four dimensions. At high 
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temperature substantial changes in the bulk parameters together with strong hadronic cor- 
relations make a simple assessment of the instanton contribution to the pressure difficult. 
The screening, however, is found to decrease at high temperature. Our estimates in matter 
becomes unreliable at temperatures of the order of T ~ m\ < N > /V3. At these tem- 
peratures, a Kosterlitz-Thouless transition in the instanton language may occur, followed 
by a substantial drop in the topological susceptibility, since x(q = 0) ~ m < tptp >~ 0. 
This transition is likely to be £7a(1) restoring and may be related to the chiral restoration 
transition in QCD. This point can be clarified by studying the static correlations in the 
topological charge at large spatial separations. 

In so far, our attitude has been to constrain the instanton effects by working from 
the long wavelength limit, where the physics is known. Alternatively, one can try to 
start from the short wavelength limit, and build up the correlations from the original 
fermion determinant in the QCD action J3J. We believe that our arguments provide 
helpful constraints in the long wavelength limit both in the vacuum and in matter. 
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